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Abstract 

I examine the arguments which have been given for quantum fluctuation-dissipation 
theorems. I distinguish between a weak form of the theorem, which is true under 
rather general conditions, and a strong form which requires a Langevin equa- 
tion for its statement. I argue that the latter has not been reliably derived in 
general. 

1 Introduction 

A quantum generalization of the Nyquist dissipation-fluctuation relation was 
discussed by a number of authors many years ago. I shall refer explicitly to the 
papers of Callen and Welton JJ and of Senitzky [2] , and to the books by Landau 
and Lifschitz [3], Feynman and Hibbs [4], Gardiner and Zoller [5], Kogan [6], 
Balescu f7| and Kreutzer [8] . I will distinguish between the weaker version (not 
including a Langevin equation) of the relation flrst discussed by Callen and 
Welton and the stronger version flrst formulated by Senitzky, and examine the 
range of validity in each case. The latter author seeks to establish a quantum 
Langevin equation containing a dissipation term and a noise term, connected 
by the fluctuation-dissipation relation. The text books [3] , [6] , [7] , follow [l] . 
The monograph ^ includes a version of [2^ . I also briefly review the formalism 
of gj. 

Another approach, which should contain much of the same information is 
based on the master equation for the density operator (see [5]). I discuss this 
approach in section 4. 

A simple quantum system is envisaged in interaction with a bath (the "loss 
mechanism") which has many degrees of freedom, closely spaced energy levels, 
and is initially in thermal equilibrium at temperature T. The bath dissipates 



energy from the simple system and at the same time feeds noise into it. The 
generahzed Nyquist relation connects the noise with the dissipation. 

There is a special case in which all the equations of motion for the bath 
degrees of freedom are linear (that is, the loss mechanism consist of oscillators 
which are linearly coupled to the simple system). In this case, the theorem is 
certainly well understood [5] . This special case covers many important applica- 
tions, such as when the bath consists of photons or phonons (if the interaction 
with the system is limited to emission and absorption). But a bath consisting 
of fermionic or spin degrees of freedom would not fall within the special case. 

The epithet "quantum" implies that terms of order h are retained in the 
expectation value of two noise operators. For consistency, this implies that the 
noise term is indeed a non-commuting operator. Since the Langevin equation 
gives the behaviour of the simple system in terms of the noise, the simple system 
must also be quantum. It has been stressed in [lOl that the expectation values 
of the noise operators in their two different orderings have different physical 
interpretations; so it is impossible to ignore the fact that they do not commute. 

The quantum Langevin equation has been quoted in connection with exper- 
iments [n] [12] (see also [6], [13]) in which a Josephson junction is shunted by 
a resistor. The Josephson phase difference 6 is modeled as the simple system, 
and the resistor is modeled as the bath. The resistor R tends to dissipate the 
Josephson junction voltage U = -^S because of its conductance 1/R. The noise 
current which couples to U must be bilinear in the electron creation and annihi- 
lation operators. Therefore the equations of motion for these electron operators 
are not linear in the complete set of degrees of freedom {U together with the 
electron operators). So it is not obvious that the Langevin equation can be 
justified. 



2 The method of Callen and Weston 

In this section, I review a weak form of the fluctuation-dissipation theorem, 
which is proved (originally in [T]) under rather general conditions, provided 
only that departures from equilibrium are small. 

The complete system has quantum variables qi,Pi with i — 1, N with N 
large. If it is driven by an external classical driving force F{t), the Hamiltonian 
is 

H{q,,pi)-F{t)Q{q,,pi), (1) 

where Q is some particular combination of the dynamical variables. (In one 
example, F is an applied voltage and Q the resulting current.) The whole 
system is, in the absence of F(t), in thermal equilibrium at temperature T. The 
eigenvalues of H are £^„, and are assumed to be closely spaced (so sums over 
n can be approximated by integrals). Thermal expectation values are given by 
the distribution function 
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Following Landau and Lifshitz [3], we expect, provided that F is small 
enough, a linear response of the form 

POO 

(Qpit)) = tiipQp) - / dt'xit - t')F{t'), (3) 
Jo 

(the suffix F on Qp indicates that this is the response to the external force F; I 
shall use Q with no suffix to denote a fluctuation in the absence of F) . Equation 
(3) defines a "linear response function" or "generalized susceptibility" x (the 
inverse of the impedance). Note that it is the expectation value oi Qp which 
appears in (3), so the whole equation refers to classical quantities. 
For the case where 

= »[Foexp(iu;t)], (4) 
we have (in terms of Fourier transformed function, denoted by " 

iLu{QF{cu))=xiLu)Fo, (5) 

where 

xi-Lu) = -n^)- (6) 

The power provided by F and dissipated into the system is 

W = (dH/dt) = -{QF)F{t) (7) 
and its mean value (over time) is 

W=i5?{x(a;)}|Fop. (8) 

But we can also calculate W by applying lowest order perturbation theory 
to the perturbing Hamiltonian —QF. The result can be expressed in terms of 
the expectation value 

S{uj,T) = ^ [ dtexp{-iujt)ti{pQ{0)Q{t)} (9) 
2n J 

Comparison of this result with (8), yields the theorem 
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huj 



exp(huj / kT) — 1 

If we choose instead in (9) the symmetrised product of operators, we get 

S,y„, = l[S{cj,T) + S{-u;,T)] = -N{lj,T)^{x{co)}, (11) 



where 

N(uj) = huj ] — ^ h - 

is the Bose distribution function. 



(12) 
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I call (10) or (11) the weak fluctuation-dissipation theorem. It is very general 
and simple, but it has the a limitation, which I shall now explain. If Q were a 
classical quantity, one could have, more generally than (3), 

/•oc 

Qf= / dt'x{t~t')F{t'). (13) 

Then it would be natural to define a fluctuation / in the force related to the 
fluctuation Q by 

/•oc 

Q{t) - / dt'x{t^t')f{t'). (14) 
It is then possible to deduce from (10) and (14) a formula for 

^f^^J diexp(-zc^Otr{/(0)/(t)}. (15) 

But in the quantum case, using only (3), one cannot take these steps, and so 
cannot go beyond the fluctuation-dissipation relation (10). 

Some authors (for example Callen and Welton 1, Landau and Lifshitz [3], 
and Kogan [6]) tacitly assume (14) even in the quantum case. Landau and 
Lifshitz write "It is convenient to write the formula [(14)] as if Q were a classical 
quantity" (my italics). These authors also call / a "fictitious random force". 
But, in my opinion, / is not fictitious; it is a real dynamical force exerted on 
the degree of freedom Q by the other degrees of freedom in qi^Pi- In order to 
exhibit this, and to justify (14) in the quantum case (if that is possible), one 
needs a Langevin equation, and this is the subject of the next section. 



3 Senitzky's derivation of a Langevin equation 

In this method, one degree of freedom Q,P is singled out, representing a simple 
system and this interacts weakly with a bath with degrees of freedom qi^pi, 
i = 1, N. It would be possible for the qi^pi variables to be quantum, while 
Q, P was a classical degree of freedom; but I am concerned with the situation in 
which Q, P are quantum as well. For quantum effects in the simple system to 
be important, it presumably must be microscopic (for instance a single atom), 
perhaps mesoscopic (like a very small Brownian particle) , or perhaps the phase 
of a superconducting condensate ([11] [T2]). 

For simplicity taking the simple system to be an oscillator, the Hamiltonian 

is 

H = ^ l^f^2g2 ^ HBiqr,p.) + H,„t{Q, P;q^,Pi) " QF{t), (16) 

with the interaction 

Hint = -aQK{qi,pi). (17) 

Here -F is a classical driving force, a is a small parameter, and K is some 
function of the bath variables. The suffix B stands for "bath". The choice Q 
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rather than P in (16) is somewhat arbitrary, but it is important that this term 
is Unear in Q, P. Some parts of the argument would go through for a general 
potential V{Q) instead of The Hamiltonian H is positive provided 

that 

Note that I have used the same notation as in section 2, but it refers to a 
slightly different physical model. 

The Heisenberg picture is used. The equations of motion are 

P = -mn^Q + aK + F, Q = —P. (18) 

m 

The special case in which Hb is quadratic in its variables and K is linear (so 
that the whole assembly is made up of harmonic oscillators and all the equations 
of motion are linear) is simple and well-studied (see for instance [5]), and there 
is no doubt about the validity of the fluctuation-dissipation theorem. But it 
is not clear that this model covers all physical examples. The argument of [5] 
however claims to be more general. 

It might be argued (see for example [4|) that, close to equilibrium, we can 
expand potentials about their equilibrium values and a quadratic approximation 
for H is necessarily justified. But if K is itself a quadratic function, then this 
approximation would decouple the simple system from the bath altogether. For 
example, if any of the bath variables are fermionic, then K cannot be a linear 
function of them. 

We define operators with a super-fix (0) to satisfy the equations of motion 
when a — 0, and to coincide with the true operators at some initial time to; so 
for example 

K'^"^{t)=K{t;a = 0), K^^^^t = h) = K{t = to). 

p(o) = -ml72g(0), mQW = P^; P^°\to) - P(to), Q^°Hto) = Q(io). (19) 

By using the Heisenberg equations of motion, Senitzky (equation (16) of [5]) 
obtains 

K{t)^K("\t) 

-T2 f dt' f dt"U{t' - t)[K{t'), [K{t"), HB{t")]Q{t")]U{t - t'), (20) 
h Jto Jto 

where 

U{t' -t)^ exp{-z(i' - t)H^°^/h]. (21) 

Neglecting terms of higher order in a, we replace K by K^^'^ and Hb by H^^^ 
on the right hand side of (20). We may also neglect the commutators of Q with 
and with since these are relatively 0{a). Using the time translation 
operator U (I differ from [2] slightly at this point), we finally get 

K{t) w i^(o)(t) ~ia [ dt"C{t,t")Q{t"), (22) 

Jto 



5 



where (with t = t + t" - t') 

C{t,t") ^^£dT[K^°\t),k^<'HT)] = ~^[K°{t),K°{t")]. (23) 

Equation (22) may now be inserted into (18) to obtain the approximate 
hnear equations for Q and P: 

P{t) = -mn^Q{t) + f{t) + aK°{t) - io? [ dt"C{t, t")Q{t"), 

J to 

Qit) = lp(t). (24) 

Here the K'^ term may be interpreted as quantum noise, and the C term contains 
a dissipative part. 

At this stage, C in (23) is a quantum operator, not a c-number, and it is 
in general a function of the two variables t, t" not just their difference. In the 
special case when Hb is quadratic and K is linear (the oscillator case), C is a 
c-number and is a function only of t ~ t" , and is independent of T. But the 
question is whether (24) can be simplified for a general Hb and K. 

At this point, Senitzky [2J makes the approximation of replacing C by its 
expectation value: 

C{t - t") « tr{pC{t, t")} = c{t - t", T), (25) 

where now 

^ exp(-g<°)/fcT) 
tTexp{-H'^°'^/kT) 

Notice that this p is not quite the same as that defined in the previous section 
in (2). 

In justifications for this approximation, Senitzky |2j writes "we ignore the 
quantum-mechanical properties of the loss-mechanism [that is, the bath] in 
terms of higher order than the second" , and "all our final results (but nor nec- 
essarily the intermediate steps) will be expectation values with respect to the 
loss- mechanism; and since the term affected involves only second and higher or- 
der interactions, only the higher order quantum mechanical effects are neglected 
in the final result" . I have not understood the distinction between higher order 
effects (higher order in a) and higher order quantum mechanical effects. It does 
not seem to me that the terms neglected in the approximation (24) are higher 
order in a or in ?i than the terms retained. 

//we make the approximation (25) and insert it into (24), we obtain an 
equation of Langevin type. If we choose to — —oo and take Fourier transforms 
(denoted by Q etc), we obtain 

[m{Q^ - uj^) + G{u}, T)]Q{uj) = iujZ{uj, T)Q{uj) = ak{uj) + F(w), (27) 
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where 

G(w,T) ^a^ [ du;'— 5(w',r), (28) 

J uj' — u} — le 

where c is the Fourier transform of c defined in (25). So, in this approximation, 
there does exist a c-number impedance Z and Z is the inverse of the suscepti- 
bihty X defined in section 2 in equation (5). The real part of Z is the resistance 
giving the dissipation. But, without the approximation (25), the C term in the 
equation of motion is a quantum operator. 

The term aK in (27) is interpreted as noise and plays the role of the force 
fluctuation / hypothesized at the end of section 2. The noise spectrum is (choos- 
ing the symmetrized product) 

2 p 

Sill{u;,T) = ^J dteM~^^t)i,{p{K{Q)K{t) + K{t)K{Q)}, (29) 

where the superfix (/) denotes the noise in / as opposed to Q. Then the 
fluctuation-dissipation theorem {if the approximation (26) were valid) would be 

S%1[uj,T)^-[^Z{lo,T)]N{T,u), (30) 

where N is the Bose function defined in (12). The theorem could alternatively 
be expressed in terms of the fluctuation Q produced by the noise by using 

iijjZ{uj)Q = aK{uj), 

where Z is defined in (27). 



4 The argument of Feynman and Hibbs 

A fluctuation-dissipation theorem is derived in section 12.9 of [1] (see also [H]). 
The argument is very similar that of Senitsky except that it uses Feynman's 
path integrals instead of Heisenberg equations of motion. It assumes that the 
Lagrangian is of second degree (so that the bath degrees of freedom can be 
explicitly integrated out). Feynman and Hibbs appear to argue that such a 
second degree Lagrangian is very general. 



5 The master equation 

A different approach to the same physical problem is via the master equation 
for the density operator p (see for example (5j chapter 5). For a critique of this 
method in the quantum case see [5]. 

Using now the interaction picture (instead of the Heisenberg picture), the 
density operator obeys the equation 

np[t) = -i[H,^t{t),p{t)], (31) 
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where Hint is defined in (17). Iterating this equation, 

p{t) = -J[i7,„t(<),p(to)] I dt'[H,nt{t), [H,^tit'),p{t')]]. (32) 

We assume the initial condition that the density operator factorizes into a bath 
part and a system part 

p{to) =^ Pb(^ PS- (33) 

We would normally choose pB to be the Boltzmann distribution. 

From (32) we can deduce equations for expectation values of operators. For 
example 

^tr{p(i)P(t)} = -mnhv{p{t)Q{t)} + atv{p{to)K{T)}Pit) 

+ '4- f dt'tr{p{t')[K{t),K{t')]Q{t')}, (34) 

^tr{p(i)Q(t)} - l-ti-{p{t)P{t)}. (35) 

//the loss mechanism consists of a set of oscillators, and if K is linear in these 
oscillator operators, then the commutator [K{t), K{t')] in (34) is a c-number. 
In that case, the last term on the right in (34) contributes a dissipative term 
to the equations of motion for the expectation values of Q and P, consistently 
with the expectation value of the Langevin equation (24) in section 3. 

In many cases (for instance for the oscillator model), the expectation value of 
the noise term (the second term on the right) in (34) vanishes. But we may look 
for the effect of noise by taking expectation values of products of two operators. 
For example, we could obtain a differential equation like (34) for 

tr{p{t)P{t)P{t)}, (36) 

and this would include expectation values of terms bilinear in K which would 
not vanish. However, from solutions of the Langevin equations (24), one could 
get (in the Heisenberg picture) expectation values like 

tl-{pHPH{tl)PHit2)} (37) 

(where the suffix H emphasizes the Heisenberg picture). Such unequal time 
correlators seem to appear less naturally in the master equation formalism (see 
for example section 5.2.1 of [5]). 

The question at issue in this paper is whether one can obtain an approximate 
Langevin-like equation when the commutator [K{t), K{t')] in (34) is not a c- 
number. This would be possible if the density matrix factorized (approximately) 
at all times: 

p{t) ^ psit) (E> PB{t), (38) 
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where the suffices refer to simple system and bath respectively. To investigate 
the validity of this approximation, we can insert (38) into (31) and trace over 
the bath variables, to get (with the Hamiltonian (16)) 

npsit)^-itrBiPBK)[Q,ps]. (39) 

In many cases, in particular in the simple case when the bath is a set of oscillators 
and K is linear in these variables (which is the simplest possible case, and for 
which the method of section 3 succeeds) the trace on the right hand of (39) is 
is zero. This means that ps{t) ~ Ps{to)- If it were legitimate to insert this 
approximation into the last term in (32) , we would get no integral equation for 
Psit) at all. Thus the approximation (38) seems to be unreasonable. 

6 Comments 

As mentioned in the Introduction, the quantum fluctuation-dissipation relation 
has been referred to in connection with experiments which might be thought 
to be sensitive to vacuum energy [llj jl2J. I have argued in this note that one 
should not appeal the quantum fluctuation-dissipation relation without checking 
if it is valid for the case in point. 

What is the interpretation of of the Bose distribution function (12) appear- 
ing in the fluctuation-dissipation relations? First of all, it is present only if the 
symmetrized product (11) is used (see [TU|). In the ideal case where the equa- 
tions of motions are linear, and so the commutator C in (23) is a c-number and 
the dissipation is independent of T, the Bose function (12) gives the complete 
temperature dependence of the noise spectrum (30). In this case, there is little 
doubt that (10) corresponds to the oscillators in the bath. If these are three 
dimensional, we would expect the dissipation in G in (27) to be proportional to 

Beyond this special case, the situation is less clear. I have argued above 
that the Langevin equation has not been proved in general. But even if we 
assume that the approximation (25) is justified, the factor N in (12) does not 
give the T-dependence of the noise spectrum, because the dissipation itself will 
be in general T-dependent. //the fluctuation-dissipation relation (11) had been 
justified generally, it would cover cases where the dissipation and noise were due 
to a bath of fermionic systems, and then the Bose factor N in (12) could not 
possibly represent the physics of the bath. 
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